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Abstract 

We explore the probability v(n, r) that a permutation sampled from 
the symmetric group of order n uniformly at random has cycles of lengths 
not exceeding r, where 1 < r < n and n —> oo. Asymptotic formulas 
valid in specified regions for the ratio n/r are obtained using the saddle 
point method combined with ideas originated in analytic number theory. 
Theorem 1 and its detailed proof are included to rectify formulas for small 
r which have been announced by a few other authors. 


1 Introduction 

The history on enumeration of decomposable structures missing large compo¬ 
nents starts from the papers by K. Dickman and N.G. de Bruijn dealing with 
natural numbers composed of small prime factors. After numerous continua¬ 
tions, this analytic theory is now extensively developed and well exposed in the 
book by G. Tenenbaum m and in more recent papers. By analogy, a similar 
theory was carried out for polynomials over a finite field (see, for example, |T6] , 
0) and generalized to the so-called additive arithmetical semigroups (see [21] , 
m , ED- The survey 0 discusses the parallelism between the theories. In no 
way, the list does not pretend to be complete, however, it has influenced the 
present paper devoted to permutations. So far, the results on this particular 
class of structures do not reach the level of research achieved for natural num¬ 
bers. We focus only on permutations comprising the symmetric group S„ and 
seek asymptotic formulas for the probability ^(n, r) that a permutation sampled 
uniformly at random has cycles of lengths not exceeding r, where 1 < r < n, 
reN, and n —> oo. The goal is to cover the whole range for the parameter r. 

Let us start from an exact formula. Denote No = N U {0}, £ r (s) = lsi + 
••• + rs r , £(s) = £ n (s), where s = ( si,...,s n ) G Nq. If kj(a) equals the 

0 AMS 2000 subject classification. Primary 05A15, secondary: 60C05, 11B75, 60F10. 
Key words and phrases. Random permutation, cycle structure, involution, saddle point 
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number of cycles in a permutation a £ S„ of length 1 < j < n and fc(cr) := 
(Aq(cr),..., k n (a)) is the cyclic structure vector, then (see, for example, j3j) 

” 1 

|{cr e S„ : k(cr) = s}| = 1 {£(s) = n}n\ —--. 

3 3 S i' 

Hence 

l _ r i 

v(n,r) =—\{a £S n : kj(a) = 0\/j £r+l,n}\= V TT —— 
nl - L - L 7 s JSi! 

i r (s)=nj=l J 3 

where the summation is over the vectors s £ Nq with £ r (s) = n. The formula 
can be rewritten in terms of independent Poisson random variables (r.vs) Zj, 
\ < j < n, given on some probability space (Q, J 7 , P) and such that E Zj = 1 /j. 
Namely, 

v{n,r) = expj^ -\p{£ r (Z) = n), (1) 

l j=l J J 

where Z := {Z \ 1 ..., Z n ). In two trivial cases, we have v(n,l) = 1/n! and 
is(n,n) = 1. It is fairly tedious to extract information from the exact formula if 
r is large. Let us discuss asymptotical behaviour as n —> oo. 

A historical overview may be started from the number of involutions in S n . 
Namely, in 1955 L. Moser and M. Wyman m proved that 

n\v(n,2) = -^=n ra/2 exp j — ^ +n 1/2 - i}(l + o(l)). ( 2 ) 

H. Wilf included a detailed proof of © into Chapter 5 of his book [22]. 
However, Exercise 8 in it gives an erroneous expression for r = 3. It shoud be 

n\v{n, 3) = ^=- exp | - y + y 2/3 + y 1/3 - y } (l + o(l)). (3) 

As we have been able to check, the last formula without a detailed proof 
firstly appears in A.N. Timashov’s paper [20]. Note that his reference to V.N. 
Sachkov’s paper m , in which formulas (23) and (24) really concern n\v(n,r ) 
for an arbitrary r, is misleading. These formulas have been also presented with¬ 
out a proof; containing a misprint, they go in contrast to © and even to ©• 
We have to note that, a year later, M. Lugo [12] also gave m leaving for a 
reader other cases of n\is(n, r). Recently E. Schmutz kindly gave a reference 
to the manuscript by T. Amdeberhan and V.H. Moll t 2] dealing with the same 
problem. Their Theorem 8.1 also contains errors. We feel obliged to present a 
correct formula in Theorem 1. 

Let r(^) be the Euler gamma-function, where z £ C. Avoiding numerous 
brackets, instead of O(-), we will use a complex quantity B 1 not the same at dif¬ 
ferent places but always bounded by an absolute constant. Otherwise, stressing 
dependence on a parameter v in an estimate, we will write O v {-) with the extra 
index. 
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Theorem 1 . If 2 < r < log n, then 


n\u{n, r) = -L n n(1 ~ 1/r) exp j ^ d rJV n (r “ JV)/r ) (l + Hn" 1/r ). 
* r jv=o ^ 

Here d r o = — 1 + 1/r, 


and 


if 1<N <r — 1. 


d r N = 


rfr,r = --V- 

r H J 


T(N + 7V/r) 

(r - 7V)r(iV + l)r(l + iV/r) 


Our main results are the next two theorems. We prefer to present them as 
asymptotic formulas for the local probability P(l r {Z) = n) having the Cauchy 
integral representation 



where a > 0 is to be chosen. In the saddle point method, we take a = x := 
x(n/r), where the function x{u ) := x r (u ) is the unique positive solution to the 
equation 

r 

x(uY = ur, u> 1. (5) 

j =i 

Evidently, the problem concerns asymptotical behaviour of the nth power series 
coefficient of the function exp j J2j< r z J Vj}> belonging to the so-called Hay- 

rnan’s class of admissible functions (see [10]). B. Harris and L. Schoenfeld [9] 
extended Hyman’s methodology in obtaining further asymptotical terms. In 
particular, it yields 

P{MZ) = n)= Jt\{x) ( X + ° r (n)) (6) 

for arbitrary bounded r. Here x = x(n/r) and 

Q( z ) : = ^ ex pjl]^—-|> \{z)~Y,jzP 

Actually, we owe to E. Schmutz whose Theorem 1 and the facts presented below 
it in |18j clarify the use of the general and fairly complicated expansion given 
in J5]. A.N. Timashov [20] mentions a Sachkov’s result from 1986, extending 
formula (0 for r = o(logn). Unfortunately, we failed to find a relevant paper. 
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The above mentioned results deal with the case when the ratio n/r is large. 
In addition, there exists a vast literature dealing with the case when n/r is 
small. In fact, the problem is related to the limit distribution of the longest 
cycle length (say, L n (a)) and other statistics of a £ S n . So, V.L. Goncharov’s 
result [B] from 1944 shows that 

v(n,n/u) = — r I{cr £ S„ : L n (cr) < n/u}\ = p(u) + o(l) 
nl 1 ' 

uniformly in u > 1. Here p{u) is the Dickman function defined as the continuous 
solution to the difference-differential equation 


up (u) + p{u — 1) = 0 

with the initial condition p(u) = 1 for 0 < u < 1. Since p(u) < T(u + 1) _1 , the 
error in the last estimate can dominate if u — > oo. Theorem 4.13 in [3], applied 
for permutations, deals with the relative error. Namely, it shows that 


v(n,r) = p{u)(\ + o(l)) 

if n/r —> u £ (0, oo). As a byproduct of enumeration of elements in an additive 
arithmetical semigroup missing large factors, the last relation (extended to a 
larger region for n/r) has appeared in the first author’s paper [14]. The result 
is contained in Theorem 3 below. The present paper fills up missing details in 
its sketchy and indirect proof. 

Theorem 2. As above, let x = x(n/r). Then 


P(e r {Z) = n ) 


Q(x) ( 1 + Br\ 

\/2ttX{x) V n ) 


provided that 1 < r < cn(logn) 1 (loglogn) 2 , where c = l/(127r 2 e) and n > 4. 

Actually, the result holds for all 1 < r < n. The bound for r is left to show 
the limitations of the applied approach. We will prove in Lemma fill of Section 
5 that 


x 


i/r _ 1 _ _ T(N + (N - 1 )/r) _ 

r N^2 ( N ~ 1 ) ]? ( 7V + 1)T((N - 1)A) 

+ I n -i+iA + 

r n 


n -( N -l)/r 


(7) 


if 1 < r < logn. Nevertheless, a direct proof of Theorem 1 by the use of © and 
Theorem 2 would be rather involved. We press more on the Lagrange-Biirmann 
inversion formulas, instead. Lemmas [4] and [5] below provide approximations 
of x and \(x) for larger r. Then it is preferable to apply another technique 
giving even a sharper remainder term. To catch an idea, let us obtain explicit 
expressions of the main term of the probability examined in Theorem 2. 
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(8) 


For u > 1, define £ = £(it) as the nonzero solution to the equation 

e 5 = 1 + uf 


and put £(1) =0. Denote also 


I(s) 




V 


s G c. 


Let 7 denote the Euler-Mascheroni constant. 


Corollary 1. If n > 4 and y/n log n < r < cro(logn) 1 (loglogn) 2 , then 




Bn\og(n/r) Br\ 
r 2 n / 


= ^(7)1- 


Bn\og(n/r ) £?r\ 

r 2 n J 


Consequently, the approximation involving the Dickman function holds in a 
rather wide region for n/r. This makes us even more greedy and motivates in 
searching for another approach to refine the corollary. 

Theorem 3. If y/n log n < r < n and n > 1, then 

Bn\og{n/r + 1)' 


P 


(£ r (Z) = n) = 1 + 


Having in mind an analogy with number theory, one can guess that the last 
approximation using Dickman’s function is hardly further extendable. The next 
corollary justifies one expression for all possible r. 


Corollary 2. For 1 < r < n, we have 


P(f r {Z) = n) 


Q(x) L + 

V / 27tA(x) \ n ) 


(9) 


On the other hand, Theorem [3] gives a better remainder term than Br/n if 
n 2 / 3 log 1 / 3 n < r < n. 

The paper is organized as follows. Section 2 collects known and new auxil¬ 
iary properties of the involved functions and the saddle point approximations. 
Theorem [2] and Corollary |T| are proved in Section 3. Section 4 is devoted to 
Theorem [3] and Corollary [2] A detailed proof of Theorem 1 is given in the last 
section. 


2 Auxiliary Lemmas 

Throughout the section, we assume that r > 2 if it is not indicated otherwise. 
Let £(w), p(u), and x(u) := x r {u) be the functions defined above for u > 1. 
Recall that they are positive and differentiable if u > 1. We will often use the 
abbreviation / = f(u) and /' = f\u) for the values at the point u, where f(v), 
v > 1 , is any of the involved functions. 
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Lemma 1. If u > 1, then logit < £ := £(it) < 2 log u, 


f — log u log log(it + 2 ) 


B loglog(u + 2) 
log(it + 2 ) 


and 


C := £'(«) 


1 g _ 1 c „ ( B ] 

u £ — 1 + 1/it it eXP \ log(it + 1) J ’ 


(10) 


Proof. To establish the effective bounds for all u > 1, it suffices to employ 
the strictly increasing function I'{v). Indeed, the lower bound follows from the 
inequality 


u = I\f{u)) 


f e^dt > P (log i 

Jo 


it — 1 

logit 


following from it log it — it + 1 > 0 if it > 1. To prove the upper estimate, it 
suffices to repeat the same argument. 

The asymptotical formulas for £(it) and its derivative can be found in fi ll 
or in the book m- 

The lemma is proved. 


Lemma 2. For it > 1, 


P(u) = yi^sxp {7-«£ + /(£)} (i + f)- 

Proof. This is Theorem 8 in Section III.5.4 of m- The result has been 
proved by K. Alladi jT]. 

Lemma 3. Let 

pOO 

p(s) := / e~ sv p(v)dv = exp {7 + /(—s)} , s £ C, 

Jo 

be the Laplace transform of p(v), s = —£(it) + ir =: —£ + it and r G R. Then 


Pis ) 


I? exp {/(£) — r 2 li/27T 2 } if |r| < 7T, 

13 exp {/(£) — u /( tt 2 + £ 2 )} if |r| > 7T 


and 


p(s) = - 1 + 


B(l+fu) 


if |t| > 1 + uf. 


Proof. This is Lemma 8.2 in Section III.5.4 of [ T9] . 

Denote a A b := min{a, b} and a V b := max{a, b} if a, b G R. Recall that 
x := x(u) is the solution to the saddle point equation and X(x) = - 

Lemma 4. If it > 3, then 


x = exp 


log ( 1 t(r A log w)) 
r 



( 11 ) 
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If 3 < u < e r , then 
x = exp 


= exp 

Moreover, for u > 1, 


| log (it log it) 

{;} 


1 + 


B log log u B log u 
r log u r 2 


1 + B log log u + B log u 


r logit 


|A(cc)/(r 2 u) — 1| < log 1 u. 


( 12 ) 


(13) 


Proof. By definition, x > 1 and u < x r < ru for u > 1. The well-known 
property of geometric and arithmetic means 


c (r+l)/2 = {x l x 2'"X r )l/r<lY,X j = 


1=1 


yields 


u 1/r <x< it 2/(r+1) < u. 

We have from the definition that 

x r = 1 + ru( 1 — aT 1 ). 

Consequently, by ta and by virtue of 1 — e t >te 1 if t > 0, 

x r > ru( 1 — exp{—(logu)/r}) > ulogitexp{—(logu)/r} > e _1 ulogu 
provided that r > logu. Similarly, 

x r < 1 + ru( 1 — exp{—2(logu)/r}) < 1 + 2ulogu. 

The last two inequalities imply 

r log x = log (u log u) + B 

for r > logu. 

If r < logu, we have 

x r > ru(l — exp{—(logu)/r}) > (l — e~ 1 )ru. 
and x r < 1 + ru. Now 

rlogx = log(ur) + B. 

The latter and (TlCJll lead to relation (fill) . 

To sharpen CD for 3 < u < e r , we iterate once more and obtain 


(14) 

(15) 


(16) 


r log x = log 
= log 


l + ru(l — x 1 ) 
1 + ru ( 1 — exp 


— log(u logu) 


(* + ?) 


log (u log(u log u) + Bu + B{u/r) log 2 u'j 
B log log u B log u 


= log(ti log u) 


logit 
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This is the first relation in m- Comparing the result and Lemma [lj we have 
the second one. 

To prove in, we first observe that 


A(a;) 


rx r+1 
x — 1 


x r+1 

(x- 


— X 

w 


7 * 4-1 

rx ^ — ru o 

---= r u + 

x — 1 


r(x — u) 
x — 1 


(17) 


Further, 

„ 1 u — x 1 1 

0 <- 7 < - 77 ^ 1 -> 

ru x — 1 r(x — 1) log u 

due to (1X41) and r(x — 1) > r(e0 og “)/’~ — 1) > logw. 

The lemma is proved. 

Using properties of differentiable functions, we improve the remainder term 
estimates. 


Lemma 5. If 1 < u < e r , £ := £(u), and := f(u), then 

X) £log(u + l) 


= x(u) = exp | -1 + 


and 


U») = Ui + 

x r V 


£'/ B log(it + 1) 


(18) 


(19) 


Proof. One may skip the trivial case when r is bounded. From © and (1131) . 
for the differentiable function x(v), we have 


„ ,, , rx(v) 

0 < X (v) = ———r < 


c(w) 


Bx{v 

\(x(v)) ~ rv(l - log -1 3)” rv 


( 20 ) 


if v > 3. The same holds if 1 < v < 3. Indeed, in this case it suffices to apply 
the trivial estimate X(x(v)) > r 2 /2 > r 2 v/ 6. 

As a function of v, exp{£(i>)/r} is also strictly increasing; therefore, given 
any u > 1 and the value £ = £(it), we can find w > 1 such that 


x(w) = exp {£/r}. 


Now 

x — exp {£/r} = x[u) — x(w) = B(u — w)x'(v), (21) 

where v is a point between the u and w, irrespective of their relative position 
on the real line. 

Using USD with w instead of it, we have 

x r (w) — 1 = e^ — 1 = nu(l — x(w)~ 1 ) = ric(l — e“^ r ). 

By the definition of £ and Lemma 1, we obtain from the last relation that 
u£ = rw{f/r + B(f/r) 2 ) with \B\ < 1/2. Hence 

\u — w\ < wf/(2r). (22) 
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If u < 3 and r > 1, then 0.09u> < w(l — £(3)) < 2u < 6 and u — w = Br 1 . 
Therefore, estimates (l20l) and (I2T1) imply 

x — exp {£/r} = Br ~ 2 , 

as desired in (ITS]) . 

If u > 3, then by virtue of £ ~ log u as u —> oo and log u < r, we obtain from 
(l22l) that | u — u>| < (3/4 )w if r is sufficiently large. Hence (4/7 )u < w < iu 
and (4/7)u < v < 4 u. By LemmalU this gives x{y) < x(4u) = B. Formula (fl8l) 
again follows from (USD and m- 

To derive approximation m of the logarithmic derivative, we use similar 
arguments. First, given u > 3, we define y > 1 such that x = e^ v ')/ r and claim 
that 

£ = £(l/)(l + B/r). (23) 

Indeed, if also u < e r , then an observation in the proof of Lemma H] gives us 
£{y) = rlogx < log(wr) < (6/5)r if r is sufficiently large. By the definitions 
and inequalities 


0 < 


1 — e~ 


-i = ^< 


t 2 /2 ^ 3 1 

t-t 2 /2 - ~2 


if 0 < t < 6/5, we further obtain 


= x x r - 1 = e gfa) - 1 j{y)/r = ( B£(y) 

r x — 1 t;(y) 1 — e~£(y)/ r r 


(24) 


with 0 < B < 3/2. Hence 15/14 < (5/14)w < y < u and also Ci v ) = B/v = B/y 
for all v e [ y , u], by LemmaU Inserting this and E41) into £ — £(?/) = (u — y)t;'(v) 
with some v &[y,u], we complete the proof of (E51) . 

Let us keep in mind the bound y > 15/14 and return to the logarithmic 
derivative. It follows from m and m that 


x 

x 


x r - 1 


"{x — 1) 


1 

ru 


Now, the idea is to rewrite the quantity in large parentheses via £(?/), then use 
inequality E?|) to approximate it by £ and £'. 

The inequality 0 < t~ l — (e 4 — l) -1 < 1 applied with t = £(y)/r gives 
(r(x — l))” 1 = 1 /£(?/) + B/r\ therefore, 


x\ viyy) r J ru 

Because of m, the first ratio inside the parentheses is l/(y^ , (j/)) which, by 
Lemma [l] satisfies an inequality 


i y log y - y +1 
vi'iy) ~ y^ogy 



> 0. 
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Now using (El and El, we obtain 




(*+?) 


x ru 1 + y£(y) — y \ r ) ru 1 + uf — u 
if 3 < u < e r . 

In the case 1 < u < 3, we have from (in 


< / 1+ -Blogu\ = ^/ + Blogu\ 
A r ) r V r ) 


Mx) = 


i=i 


>t/r . 


B 


= r - e U/ r _|_ B r = r 2 I te^dt + Br 


3 =1 


Hence 


? ,2 r 2 
— (u£ + 1 - u) + Br = — + Br. 


*(' + $)■ 


x X(x) r 

The lemma is proved. 

We will need an estimate of the following function 


T(z) := 


e* — 1 ft e‘/ r 


— 1 ) dt, z € C. 


t \r e*/ 1 " — 1 

Lemma 6. If z = rj + it , 0 < rj < nr and —nr < r < nr, then 


T(z) 


2 r 


< 


4e v 


12 r 2 ' 


Proof. The well known theory of Bernoulli numbers {&„}, n > 0, gives us 
the series 




(-l) fc+1 C(2/c) 


,2k 


(2tt) 


2fc 


(26) 


n —0 fc=l 

-.—2k /• /-/o\ _ _2 , 


converging for z G C. Here £(2fc) = X^ m >i m 2k — C(2) = 7r 2 /6. Hence 


oo 


rW = (e t -l)dt + 2j2 


(-l) fc+1 C(2 fc) 


2r 

_(e'- z -l) + 2 ^ 

OO 


fc=l 


(27rr) 


2/c 


[ (e* — 

Jo 


(-l) fc+1 C(2fc) 


- (27rr) 

(-l) fc C(2fc)0 2fc 


2fc 


eV*" 1 


-(2k-l)J e t t 2k ~ 2 dt^j 


k =1 


2/c(27rr) 2fc 


(27) 


Under assumed conditions, |*| 2 < 2-7r 2 r 2 ; therefore, summing up the series, 
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we obtain 


e' 1 2t r 2 Izl 2 ^ 1 7T 2 




— v 


12 k 


< 


3 (27rr) 2fe ' 6 k{2nr) 2k 

on e’Ty + |t|) 7y 2 + r 2 e v ( 2 tt i r\ 

-1-12-1_H _l J- < — M j-1-j_ 

r 3r 2 12r 2 “ r V 3 12/ 


12r 2 


Ae* 1 t 2 

— r + I2 r 2 


The lemma is proved. 


3 Proof of Theorem 2 

The non-standard part of our proof concerns the following trigonometric sum 

grit, y ) := 1 -7-) t e (~ n ’ 7 r l> y > !• 

j<r 

Its behavior outside a vicinity of the point t = 0 is rather complicated; therefore, 
we consider it in a separate lemma. Denote 

r 

A k := y^' fc -V, k> 1, 

3=1 

where x = x{n/r). In particular, Ai = ur and A 2 = A(x). 

Lemma 7. If t € [— 7r, 7r] and y > l, then 


®-gr{t,y) < - 


2 y 


r+1 


+ 


2y 


7 r 2 r(y — 1 ) [y — l) 2 + f 2 r(y — 1)' 
If 1/r < |f| < 7 t, x = x(u), and u := n/r > 3, t/ien 

1 u l-4/(r+l) 2 2 


5ft3 r (<) := 5fty r (f,x) < - 


47t 2 log 2 


logn 


(28) 


(29) 


Proof. Observe that 




3=1 


3 = 1 

y r+l ^ e »«(r+D( y -l) _ 


r(y - 1) 
„r+l 


< 


y 


r (y - !) Vlye** - 1 


ye lt — 1 
2 /- 1 


r(y -1) 


1 — 5ft 


e u {y ~ 1) 
ye 1 * — 1 


-1 + 


2 y 


r{y - 1)' 


(30) 
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If |t| < 7 r, then 


I ye u - 1| = f, 2t/(l — cost) \ 2 ((y ~ l) 2 + (4/7r 2 )f 2 )2 

I /-1 V (y- 1) 2 ) - y- 1 


because of 
Using also 


a 


2< 2 /tt 2 < 1 — cost < f 2 /2. 

1 i> 2 


Va 2 +v 2 


_ 1 < _: 


2 a 2 + u 2 ’ 

with a = y — 1 and t; = ( 2 / 7 r)t, we obtain 


a > 0, d£R, 


y -1 

|ye lt - 1| 


-1 < - 


2 t 2 

^ (y — l ) 2 + t 2 ' 


(31) 


Inserting this into (131)1) , we complete the proof of inequality 

If y = x, 1/r < |f| < 7T and u > 3, we combine (E51l with estimate (fill) . We 
have 


x r+1 x 

= n + 


x — l 


x — 1 


> ur 


and 


1 < logit < r(x — 1) < r(u 2 ^ r+1 ) — 1) <- u 2 /(r+i) logit. 

r + 1 


So, we obtain 


5R 9r{t ) < 

< 


1 u 2 / 1 \ 

7T 2 r 2 (x — l) 2 r V cc — 1/ 

/ r + 1 \ 2 it 1-4 /( r+1 ) 2 2 

V 27T7- ) log 2 it r logit' 


Lemma [7] is proved. 

Proof of Theorem [2j As it has been mentioned in the Introduction, it suffices 
to examine the case when r > 4 and n is large. In the introduced notation, we 
have u > c -1 (logn)(loglogit) 2 and 


P(t r {Z) = n ) 


Q{x ) 

27T 



^(Ri(n) + R 2 (fi)) 


exp {g r (t)} e ltn dt 


(32) 


with to := r~ 1 / 12 n -5 / 12 . 

Expanding the integrand in K\(n ), we use relations e lt = 1 + it — t 2 /2 — 
it 3 /6+Bt 4 if t G R and e w = 1 + R|it;|el u 'l if it; G C. Consequently, checking that 
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< ( r 3 n)(r 7 / 3 n 5 / 3 ) = (r/n) 2 / 3 < 1 and using the abbreviation A := A2, 
we obtain 

exp {g r {t)} = exp {iAif — (A/2)t 2 — i(X 3 /6)t 3 + -BA-jt 4 } 

= exp {iiAi — (A/2)t 2 } (l — i(X 3 /6)t 3 + BX 3 t 6 ) + -BA4C 1 exp { — (A/2)f 2 } 
= exp {itAi — (A/2)f 2 } (l — i(X 3 /6)t 3 ) + B^X^t 4 + Agf 6 ) exp { — (A/2)f 2 }. 

Recall that u = n/r , Ai = n , A& < r k u if k > 1, and, by LemmalU A = A(cc) ~ nr 
as n —> 00 because of u —> 00. We now see that 


Ki{n) 



-(A/2)t a dt + 


B 

71 




Vx, 


\v\>t 0 VX 


2 / 2 dv + 


B 

u\J A 



B 

uyfX 


Considering K 2 (n), we first observe that, by virtue of (l3ll) . 5ftgy (f) < — (2/-7r 2 )At 2 
if to < |t| < 1 jr. Therefore, the contribution of the integral over this interval 
to K 2 (n ) equals B/uy/X. 

Further, we apply Lemmas 0] and [7] to get 

Ki(n) = B max I exp I H-= 

i/r<\t[<n 1 PiyrW/l U y/X 

B [ 1 u i-4/C+i) 1 I B 

= ci exp \ _ i?n^r + 2 los " + log ''} + wi- 

It remains to prove that the quantity in the large curly braces does not exceed 
— log u + B if the bounds of r are as in Theorem [2] This is trivial, if 4 log it > 
r + 1 > 5. If 4 log it < r + 1 and n is sufficiently large, we have an estimate 

1 ■u 1_4 /( r + 1 ) > 3 cu > 31ogn(loglogn) 2 ^ 

47t 2 log 2 u ~ log 2 it “ (log log n + 2 log log log n + B) 2 


which assures the desired bound K 2 ( n ) = B/uy/\. 

Inserting the estimates of Kj(n), j = 1,2, into (l32l) . we finish the proof of 
Theorem 0] 

Proof of Corollary |T[ In the above notation, we can rewrite 


log<5(tc) = —nlogx+ / 1 dt=—n log 

Jl 3 =1 

l-rlogx v _ I 

= —nloga;+ / 

Jo 


t r - 1 
t- 1 


dt 


dv 


v r 1 — e~ v / r 
= —ur log x + I(r log x) + T(r log x). 


(33) 


If 1 < u < yjn/ logn, then itlog(u + 1 ) = Br and, by Lemma 01 
B log(it + 1) 


r log x = £ + 


I(r logx) — 4(C) = -Bit(log(it + l))/r. 
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Thus, by Lemma 0 we obtain 

T(r logo;) < 4 x r /r = Be^/r = B(u\og{u + 1 ))/r. 
Inserting these estimates into drill , we deduce 

Q{x) = exp{—u£+ !(£)} (l + Bulog ( u + 1 ^ . 

Observe also that, by Lemma 0 a relation 

1 _ [x r _VCf 1 | -Blog(u+l) \ 

V^A \ rx r \ r J 


(34) 


(35) 


holds if 1 < u < \Jn/ logn. 

Now, it suffices to apply the last two relations only for c -1 (log n) (log log n ) 2 < 
u < a Jnl logn. Taking into account Lemma El we can present the formula in 
Theorem E] in two ways: 

P(£ r {Z)=n) = -^=exp{-< +J(g)}(l+ Bulogu + 

r\/2n V r uJ 

e -7 . ,/ Hu logn H\ 

= -P(«)(l +-—+ -)■ 

r V r u / 

Corollary [l] is proved. 


4 Proof of Theorem 3 


The idea is to use the Cauchy integral 0 with a = y := e^ 1 " which is a 
good approximation of the saddle point. Here, as above, £ = £(u) is defined 
by the relation e^ = 1 + for u > 1 and £(1) = 0. Such a choice relates 
Q[z ) with the Laplace transform of Dickman’s function. Namely, if z = e~ s ^ r , 
s = —£ T irt =: —£ + it, and |i| < n, then, as in (1551) . 

Q (e _s / r ) = exp {us + I(—s) + T{— s)} = p(s) exp { — 7 + us + T(—s)}, (36) 

where T(—s) is the function examined in Lemma[ 6 ] 

Observe that, under the conditions of Theorem 2, 1 < u < ^n / logn, where 
n may be considered large. Let us introduce the following vertical line segments 
in the complex plane: 


A 0 :={s = 
A 2 := {s = - 
and Aoo = {s = 

0 

P(£ r (Z) = n ) 


—£ + ir : |r| < 7 r}, Ai := {s = —£ + ir : 7 r < r < r7r}, 

-£ + ir : — 7rr < r < — 7r}, A = {s = — £ + ir : |r| < r’Tr}, 

= — £ + ir : |r| > ?' 7 r}. Taking into account (ED. we have from 



Q{z)dz 

z 


= / A ■»”«*)*+ |A - 4* 

=: / + J. 
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Using Lemmas [I] [2l and [3] for the case |r| > nr > 1 + zt£, we obtain 
e _7 p(u) 1 


/ = 


— / /K*Me“ s ) 

iru J 


r 2iriru 

e -7 p(n) i Be~ u ^ + 1 


r wr 2 ' 2mur J A 

e~' 1 p(u) Be^~ u ^ 

~2 

r 

e~ ry p(u) Bp(u)e^~ I( -^ 
r + r* = 
e_7 p(w) / , S' 


e us p(s) e “ l ds 


■0 + 7 )- 


In the last step, we have used the fact that /(£) ~ e^/£ as £ —>• 00. 

The next task is to estimate J. If s £ A then, by Lemma [6l T(—s) = 13 
and exp{T(—s)} = 1 + BT(—s). Let us split J into the sum of three integrals 
Jk over the strips A*,, where k = 0,1 and 2, respectively. If s G A 0 then 
T(—s) = 13(1 + wlogit)/r\ Therefore, using Lemmas [TJ [H and 0 now for the 
case |t| < 7r, we derive 


Jo = 


13(1 + ulogu) 


|/5(s)e us | \ds\ 


'An 


13(1 + ulogu)p(u)y/u 


13(1 + u log u)p(u) 


*/^) d7 


Further, 


Ji — 


—— / p(s)(e T( — l)de 
Tiur 7 Ai 


2niur J A 
Be~< 


/5(—£ + 7T z)T(£ — 7ri) + 




p(—£ + nri)T (£ — 7rri)| 


3 —-J e us (^p(s)'(e T ( ^ — l) — p(s)T'(— s)e T( - S ^ds 
B 

—• L\ ~\~ L 2 H-Z/3. 

ur 

To estimate Li, we combine the first estimate of p(s) given in Lemma[3]with 
Lemmas Q] and 0 So we obtain 

= B{\ +nlogn) £+J(£) = Bp(u)(\ + u\ogu) 

1 ur 2 r 2 

Similarly, the second estimate in Lemma [3] leads to 

Be~ u ^ Bp(u)e~ 1 ^ Bp(u) 


L 2 — 


r 2 yju 
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Estimation of the integral L 3 is more subtle. It uses an estimate 
1 - b(-s/r ) - T(—s) = b(—+ - 2 ) 

\r r J 

following from Lemma [6] and the asymptotic formula b(v ) = 1 + v/2 + Bv 2 for 
Id < n\/2. We have 


La — 


ds 


j" e us - -—p(s) ^1 — b(— s/r)e T ^ _s ^ ds 

£ e-2r_zlp( s) ^(1 _ fo (- s /r) - T(—s)) + + T(-sf )) 


= Be~< 


Ie- s -l| 


' Ai 


I/K s )|(7 + 7 t)H- 


Using the two different estimates of p(s) on the line segments An := {s 6 Ai : 
|3s| < 1 + ri£} and A 12 := A x \ A n given by Lemma[H we proceed as follows: 


L 3 = B exp l —«£ + /(£) — 


7T 2 + 


+e 


1 /e ? 


/An |s| V r 


\ds\ 


+Lexp{-< + £} [ -th{— + 

Va 12 M 2 V r 


\ds\ 


D j ^ , T/ ^ u , n ,|l + ? , Be xp{-< + ^} 
= Lexp<j -n$ + /(0~ 2 ,. 2 + 2 C 


7T 2 + £ 2 


Bp(u)y/u log(u + 2) 


Collecting the obtained estimates, we obtain 

T _ r , r , 5 r Lp(w) (1 + U log u) 

— 1j\ -r L/2 - ^3 — - 9 -• 

The same holds for integral J 2 . Consequently, 

P(l r (Z) = n) = I + J„ + J, + h = AW(l + g(1 + " loe ") 

\ 7 r V r 

Theorem 3 is proved. 

Proof of Corollary [2] Combine Theorems 2 and 3 with relations (l34l) and 
(13511 valid in the region which is not covered by Theorem 2. 


5 Proof of Theorem 1 

Most of the lemmata of this section are well known and could be found in the 
literature. Let C[[u]] be the set of formal power series over the field C and let 
[u n ]g(u) denote the nth coefficient of g(u) £ C[[tt]] where n £ No. 


ds 
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Lemma 8. Let 


OO 

$(u) = ^ $nU N 
N=0 

be a power series in C[[u]] with $0 = 1. Then, the equation u = z<P(u) admits 
a unique solution 

OO 

u = f{z) = fnz N , f N = N> 1. 

N =1 

Proof. This is Lagrange-Biirmann Inversion Formula, presented, for instance 
on page 732 of a fairly concise book |4: . 

In this note, superpositions of series involving f(z ) are used, therefore we 
recall more variants of the inversion formula. Let us stress that, by Lemma [HI 
/i = 1; therefore, z/f(z) and log (z/f(z)) have formal power series expansions. 

Lemma 9. Let f(z) be as in Lemma [8] and j £ N. Then 


if N £ N 0 \ {j} and 



Moreover, 

[^o Sl ^ = ~[u N Mu) N 

if N > 1. 

Proof. Without a proof the first part of Lemma O is exposed as A. 11 on 
pages 732-733 of [33; an inaccuracy is left in the case N = j , however. For 
readers convenience, we provide a sketch of a proof. 

Let N £ N 0 be fixed. The coefficients under consideration have expressions 
in terms of $fc with 0 < k < N only; therefore, we may assume that $(w) is a 
polynomial of degree N. Then f(z) is well defined as an analytic function in a 
vicinity of the zero point. Thus, we may apply Cauchy’s formula. Afterwards 
let S and <5i be sufficiently small positive constants. Using a substitution 2 = 
u/$(u) and properties of the one-to-one conformal mapping of the vicinities of 
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the zero points in the 2- and u-complex planes, we obtain 



1 f dz 

2th J\ z \=$ f(zyz N+1 -i 


1 f d(u/$(u)) 

27T* J\u\=Si vd (u/tf>(u)) N+1 J 

1 r $('u) iv_J du 1 f <fr(u) N ~ : '- 1 d<5>(u) 

2t ri J\ u \= Sl u N+1 2m J N=(5l u N 


[u N ]$(u) N ~ j 

[u N ]^(u) N ~ j 


1 r d<I>(u) N -J 

2n(N - j)i J H=Sl u N 

N f <S>(u) N - j du 

2n{N - j)i u N+1 


j 

j-N 


[m w ]$(u) jv " j 


provided that N ^ j. 

The same argument gives 



Finally, applying the previous substitution, we derive 


[- jY ] log 7TT = b^v 7 / 

/(-) 2tt7Vi f{z) 

1 r d$( u ) N 
2nN 2 i J\ u \= Sl u N 

= -[/]$«■ 

The lemma is proved. 

We will apply the lemmas in a very particular case. Then the first power 
series coefficients of implicitly defined functions attain a simple form. Let 1{-} 
stand for the indicator function. 


Lemma 10. Let k,r,j £ N, y = y{z) satisfy an equation 



and let g{z) := z/y(z), then the following assertions hold. 
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(/) If g(z) j =: EaLo 9n ] z N , then 


9^ = T J 


E 


•Z + rj 
l , m G Nq 


? - TV ^ I l 

J rl+m = N \ 

for TV S N 0 \ {j} and 


({N - j)/ir\ i(m-l + {N - j)/r 


m 


gj >) = l{j = O(modr)} - 


(//) If log g(z) =: E“= 1 > ^en 


‘«-s £ 


rZ + m=7V 
Z,meN 0 


N/r\, .,/to— 1 + TV/r 


N> 1. 


(Ill) If 


then h_ r = 1 /r, 


and 


AW:= ?Sw? = ^ 


i=i 


N=—r 


^o = --E- 

r z —' 7 

7=2 J 


TV + r 

ft TV = —^—Ojv+r 


for TV = —r + 1, — r + 2,... and TV yT 0. 

(rtO If 


A « : = KEvw =E^* 


s,M J 


— 1 oo 


iV=0 


i/ien Ao = 1/r and Ajv = —Nb^/r for TV = 1, 2,.... 

Proof. To prove (/), combine Lemmas [8] and [5] with an equality 




y 


rl+m=N 

I,»>6No 


to — 1 + a 


TO 


, TV G N 0 , a e R. 


For (1371) . apply the second part of Lemma [9] to obtain 




= 1 {j = 0(modr)} — 


Similarly, (II) follows from the last formula in Lemma [9] 


(37) 
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Having in mind that z r h(z) has a power series expansion in C[[z]], we may 
apply the same principles. Using (l37l) . it is easy to check that 


ho — 


2ni 


f h(z)dz 11 f gizVdz 

L.— 


1=1 J 1=2 J 


Further, we observe that 


i -y r 


Hence 


y~ = 

y r (i-y) ^ 


\ v', ^ y'< \ i \ i i 

h (z) = --(z) y — = --(«)— = - W— - -?+T, 


(38) 


y 


z± 0 . 


This implies that 
= 


1=1 


1 


dh(z 


2n Ni 
1 


rlV 


M=<5 2 
f d(\ogg(z)) 


dz 


2nNi J | z!=5 z N+r 

N + r 
~^b N+r 


2irNi J | 2 | =5 z N+r + 1 


if TV > — r + 1 and h- r = 1/r. 

To prove (TV), we use relation (l38l) again. Differentiating it, we arrive at 

A(z) = = 

r y r r 


/ 00 


JV=1 

The assertion (IV) now is evident. 

The lemma is proved. 

Corollary 3. As above, let g(z) = z/y(z). Then g 0 = 1, g 1 = —1/r, 

T(N + (N - l)/r) 


9n = 


if 2 < N < r — 1, and 


9r — 


(l-N)T(N + l)T((N-l)/r) 


T(r + 1 — 1/r) 


(1 — r)T(r + l)r(l — 1/r) r 
Moreover, \gw\ < —Uj-r (JV_1 )/ r ' if N >2. 
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Proof. Apply (I) of Lemma flOl for j = 1. If 2 < TV < r — 1, the relevant 
sum has the only nonzero summand corresponding to the pair ( l,m ) = (0,7V). 
A formula for g r has two summands giving the expression. If TV > 2, then by 
Lemma [TO] and Cauchy’s inequality, 


\9n\ = 


1 


TV- 1 


[y N ](i + y + ---+y r - 1 ) (N - 1)/l 


< 


AN-l)/r 


TV- 1 


The corollary is proved. 

Corollary 4. We have 


b N = — 


T(7V + TV/r) 
TVr(TV + l)r(TV/r) 


(39) 


*/lATV<r — 1 and b r = 0. 

Moreover, 

7V|6/v| < 1 if TV < r — 1, 

= BN/r if r < TV < 2r — 1, 

and 

7V|6jv| < r N/r if TV > 1. 

Proof Again, ifl<TV<r—1, it suffices to observe that the relevant sum 
(II) of Lemma [TO] has the only nonzero summand corresponding to the pair 
( l,m ) = (0, TV). A formula for b r has two subtracting summands. 

If TV < r — 1, the given estimate follows from (1391) . If r < TV < 2r — 1, 
assertion (II) in Lemma [TUI gives 


6jv = — 


1 (N — 1 + TV/r 
N 

N 


N 

1 

r 


1/TV — r — 1 + TV/r 
TV — r 


1 + 


TV/r - 1 


n( 

k=2 




TV 

r-2 


N—r 

n( 

fc=2 


1 + 


TV/r - 1' 


— exp 


r 

/TV \ , r ) BTV 

( 7 -l)k.gAT =—. 


k=2 


We have applied an inequality log(l + x)<xifa;>0. 
Finally, by Cauchy’s inequality, 


7VM = 


[y 


N 1 


l~y r 

1-2/ 


iV/r 


= |[2/ W ](l + 2/ + ---+2/ r “ 1 ) JV/ 1<r 


iV/ r 


(40) 


if TV > 1. 

The corollary is proved. 

We now prove the promised expansion 0. 
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Lemma 11. If 2 < r < logn, then 


1 /r 

x = n 1 




T(N + (N — l)/r) 


r ^r 2 (^-i)r(iv + i)r((^-i)/r) 
_I n -i+iA + £. 


(jv-l)/r 


r n 

Proof. The equation defining a; can be rewritten as 


x 1 = 


1 — x 
1 — x 


— r \ 1 / r 


-1 


-1/r 


This gives the relation y(n~ 1 ^ r ) = a; -1 , where y = j/(z) has been explored in 
Lemma 1101 Consequently, we may apply the expansions of g(z ) given in (!) 
with respect to powers of z = n -1 / r . The first coefficients have been calculated 
in Corollary [3] It remains to estimate the remainder. Using also the obtained 
estimates, we have 


E E V 1 ,r A N < 


|r+l 


N=r+1 


N=r+1 


1 - v^e” 1 


if \z\ < e 1 . Consequently, we obtain 


1/r 

x = n 1 


E 9 Nn 


-N/r 


JV=0 


1/r L 

= nr' - 


B 


T(N+{N-l)/r) 


= _v - 

r .4-1 (AT - 1)T(N + 1)T((1V - 1) /r) 


-(N-l)/r 


N=2 


+ -n~ 1+1/r + 


as desired. 


□ 


Proof of Theorem 1. Let us preserve the notation introduced in Lemma [TUI 
First of all we seek a simple expression containing the first terms in an expansion 
of 

r i ^ 

K(z) := E ■ / w -rclog~E = H-) ~ nlogg(z). 

ffr[3y{zy y{z) 

Let D(x) := exp |X^j=i y|, we have 

log .D(a;) — n log x = A'( n_1//r ) — U '°° n ■ (41) 

r 

Define the functions R(z) and K r (z) by 

o r -1 

K{z)= E h N z N -nJ2 b N z N + R(z) = K r (z) + R(z) 

N=—r +1 N =1 
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We claim that R(z) = B{\z\ + n\z\ r+1 ) if \z\ < e 1 implying 

R(n~ 1/r ) = Bn~ 1/r . (42) 


for r < logn. Indeed, by (III) of Lemma [TUI and the estimates in Corollary [H 
we have 


EIMM W 

N=1 


r— 1 oo 

E + E 

N=1 N=r 
r— 1 AT , 


N + r 
N 


\b N +r\\z 


N 








N=1 


N—r 


B\z\ 


if 12; | < e x . Similarly, 

OO 

E \b N \\z\ N = B\z\^ 

N=r-\-l 

if \z\ < e -1 . The last two estimates yield our claim and 63. 
For the main term, we obtain from Lemma 1 101 that 


i4(n“ 1/r ) = h 0 + h Nfi~ N/r - b N n ( - r - N)/r + h_ r 

N=-r -\-1 


N=1 


r— 1 


r — N 


= h 0 -J2 ——b r - N n N/r - J2 b N n {r - N)/r + h_, 


N=1 
r —1 


n 


N=1 


= —b r _ N n N/r + h_ r 


N =1 


r— 1 


1 1 

--E-+ r E 


1 T(N + N/r ) 

--tjA 


r ^ j ^ iV(r - AT) T(IV + l)T(N/r) 


It remains to approximate 


(att) 1/2 = A a( "' 1, ’' )1/2 = T= f 1 - E Nb ""- N,r 
\\{x)/ vn \AirV ^ 


Af)/r 


1/2 


n 

r 


By virtue of Corollary[H IV|&jv| < 1 if iV < r and iV|6w| < r Ar / r if iV > 1. 
Thus, if 2 < r < logn, then 


^2 N\b N \n~ N / r < (5/2 )n~ 1/r < (5/2)e“ 1 < 1. 
N =1 


Consequently, 




(l + Bn~ x l r ). 
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We now return to probabilities. Applying ([]]), ©, (HT1) . (H2J) , Stirling’s 
formula and the last estimate, we obtain 

n\v(n,r) = ? ' ! n~ n ! r exp {K r {n~ 1/r )} (l + Bn~ 1/r ) 

\J2-n\{x ) 

ra(l-l/r) 

= -—exp { — n + Jv r .(n -1 / 7 ’)} (l + i?n _1 / r ) 

for all 2 < r < log n. 

Theorem 1 is proved. 

Concluding Remark. The approach can be adopted for more general de¬ 
composable structures, in particular, for the so-called logarithmic classes of set 
constructions (see 0)- Using a different method, X. Gourdon [7] has estab¬ 
lished some results related to our Theorem 2. They concern the asymptotic 
distribution of size of the largest component. 
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